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Abstract 

We analyse recent data on the production of forward neutrons in deep inelas- 
tic scattering at HERA in the framework of a perturbative QCD description for 
semi-inclusive processes, which includes fracture functions. Using a model esti- 
mate for the non-perturbative piece of the fragmentation process, in fairly good 
agreement with the available data, we analyse the Q 2 dependence of the result- 
ing fracture functions, which is driven by non-homogeneous evolution equations. 
We also propose a measurement of the pion production cross section in the tar- 
get fragmentation region as a new test of perturbative QCD, which in this case 
predicts also a different Q 2 evolution with respect to the one given by the usual 
Altarelli-Parisi equations. 
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1. Introduction 



During the last three decades, deep inelastic scattering (DIS) was shown 
to be the most efficient experiment designed to extract information about the 
hadron structure. With the advent of new and more powerful accelerators, it 
has been possible to make new tests of QCD in a hitherto forbidden kinemat- 
ical region (x < 1 and Q 2 3> M 2 ), obtaining very interesting results on the 
behaviour of the structure function F2, in perfect agreement with the energy 
scale evolution predicted by the next to leading order (NLO) Altarelli-Parisi 
(AP) U equations and its aproximation in the double asymptotic limit 0. 

Moreover, with the analyses of the longitudinal structure function Fl and 
the above-mentioned scaling violations, it has been possible to constrain, if 
not to measure, the gluon distribution over a large kinematical range and 
also to obtain reliable results for the strong coupling constant a s . The use 
of evolution equations has not been limited to unpolarized scattering and its 
role has been preponderant in understanding the gluonic contribution to the 
Ellis- Jaffe spin-dependent sum rules || through the U(l) anomaly Q and, 
therefore, in clarifying the situation with respect to the so-called "spin crisis" . 
In that way, a large program has been developed around the measurement 
and the theoretical analysis of DIS structure functions with successful results 
on both sides, making DIS experiments one of the most important tests of 
QCD. 

Also within DIS experiments, other interesting processes are now under 
study at HERA, specifically events of a diffractive nature or, more precisely, 
with large rapidity gaps |J. For them, the proton remains almost intact after 
the interaction with the photon, changing only slightly its momentum and 
leaving the interaction region in the very forward direction without being 
detected. This particular kind of process has been recently tested by the 
ZEUS Collaboration [R] , which has included a Leading Proton Spectometer in 
order to detect hadrons produced in the forward direction, and has measured 
the dependence of the diffractive cross section on the transferred momentum 
(t) between the initial and the final proton. In this, they have obtained results 
compatible with previous determinations, all of them in good agreement with 
the effective pomeron exchange picture, where the cross section is factorized 
into the product of a flux of pomerons in the proton (oc a;~ a ) times the, 
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expectedly universal, distribution of partons in the pomeron. 

More recently, the ZEUS Collaboration M has measured events where 
neutrons, instead of protons, are produced in the forward direction, obtain- 
ing a sizeable contribution of leading neutrons to the DIS cross section (as 
much as 10%). This opens a new window to study hard processes in a new 
kinematical region hitherto considered only in soft analyses. 

Along these experimental improvements, a new theoretical idea has been 
developed in order to deal with such kind of forward processes. In a more 
general perturbative picture of semi-inclusive processes, it has been shown 
U that both the current and the target fragmentation regions have to be 
considered. The first one by means of a cross section that can be factorized 
into a hard cross section, parton distributions and fragmentation functions. 
In the second one, with another perturbatively computable hard cross section 
but introducing new distributions called fracture functions (measuring the 
probabilities of finding a parton and a hadron in the target). 

Moreover, it has also been shown that, at NLO, fracture functions are 
essential in order to factorize collinear singularities related to the emission of 
partons in the target direction, in processes where the polarization is either 
neglected or taken into account [TT| . 



Fracture functions are expected to give the dominant contribution to cross 
sections for the production of leading hadrons in the target fragmentation 
region. In a given kinematical range, these distributions can be related to 
the parton distributions of the object exchanged between the initial and the 
final state (a pomeron in ep — > epX, or a 7r + in ep — > enX, for example); 
therefore, their measurement can be used in order to obtain structure func- 
tions of "virtual targets" , which are beyond the possibilities in a "real target" 
experiment. It is worth noting that fracture functions allow the most appro- 
priate way to factorize the hadronization properties of the forward direction 
into non-perturbative distributions and treat them in the context of QCD, 
whereas other factorization approaches represent effective pictures only valid 
in a restricted kinematical region. 

In that sense, since fracture functions obey evolution equations that differ 
from the Altarelli-Parisi ones by an inhomogeneous term proportional to 
the well-known structure and fragmentation functions, their measurement 
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can be taken as a test of new aspects of QCD, and at the same time can 
also probe the range of applicability of the cruder factorization approaches. 
Finally, since both the experimental and theoretical branches of knowledge 
have been improved during the last years, it seems that a program similar 
to one developed for structure functions can be applied to fracture functions 
with, hopefully, a similar success. 

Aiming in that direction, we analyse here the most recent experimental 
results for the production of neutrons in the forward direction in DIS, in the 
framework of fracture functions, modelling the input fracture distributions 



exploiting pion exchange ideas, as it was proposed in ref. fl2fl , and extend 
that analysis to the production of pions in the target fragmentation region. 
This last process is expected to show a clearly different evolution behaviour 
with respect to that driven by the usual evolution equations. The measure- 
ment could also be useful in order to test the non-perturbative flux proposed 
for pions and neutrons in protons. 

This paper is organized as follows, in the next section we briefly out- 



line the main features of fracture functions and use the model of ref. [12 
translated into the fracture function language for the description of the frag- 
mentation process, to analyse the data. In the third section we propose the 
measurement of forward pions and we estimate the contribution from frac- 
ture functions and from the background composed of very forward current 
fragmentation events. We also analyse there the evolution of these fracture 
functions and discuss the possibility of finding deviations from the usual 
evolution. In the last section we summarize our results and present our con- 
clusions. 



2. Perturbative treatment of target fragmentation in 
the very forward region, model estimates and their com- 
parison with the available data. 

In the quark-parton model, the semi-inclusive cross section for the pro- 
duction of a hadron h from the deep inelastic scattering of charged leptons 
carrying momentum I off nucleons of momentum P, is usually described in 



3 



terms of the variables 



Q 2 _P-q P-h_ E h l-cosg fe 

X ~2P-q' V P-V Zh p-q E p (l-x) 2 ' [) 

where q is the transferred momentum (— q 2 = Q 2 ) and Eh, E p and 9^ are 
the produced hadron and target nucleon energies, and the angle between the 
hadron and the target in the centre of mass of the virtual photon-proton 
system, respectively. Then, the most naive expression for the unpolarized 
cross section is 

do} (1 + (1 - yf 



dx dy dzh 2y 2 



°i fi/p( X ) D h/i(z h ) , (2) 



= <2:<? 



where q = 4ne 2 a 2 /x(P + I) 2 , fi/ p is the parton distribution of flavour i 
and Dh/i is the fragmentation function of a hadron h from a parton i. Next 
to leading order corrections to this cross section are also known, and have 
been shown to give a very good description of data with 9h > tt/2 JL3], the 



so-called current fragmentation region. However, the target fragmentation 
region, which corresponds to 9h = (z^ = 0), cannot be described with 
the simplified picture of eq. (2). First of all, it is easy to see that, at the 
lowest order, hadrons can only be produced antiparallel to the target nucleon 
(9h = 7r), excluding the forward configurations. On the other hand, going 
to next to leading order, the corrections to the cross section (more precisely, 
the next to leading order coefficients) develop terms proportional to l/z^. 
This divergence is not only related to soft emission (Eh = 0), but also to 
collinear configurations where hadrons are produced in the direction of the 
remnant target (9h = 0). In this way, since at lowest order hadrons cannot be 
produced in that direction, it is not possible to factorize the divergence, with 
the usual procedure, into parton distributions and fragmentation functions. 

Then, in order to describe hadrons produced in the target fragmentation 
region at the lowest order, and also to be able to perform at higher orders a 
consistent factorization of divergences originated in the current fragmentation 
region (when a parton is emitted collinearly with the target), a new distribu- 
tion has to be introduced, the so-called fracture functions, Mi t h/N{x, (1— x)z) 
[§, [|. These distributions represent the probability of finding a parton of 
flavour i and a hadron h in the target N. It has also been shown that, within 
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this picture, it is convenient to introduce a new variable z = Eh/E p (l — x), 
which allows a discrimination between the Oh = and the Eh = configura- 
tions, both leading to divergences but of a different nature. The variable z 
is equal to zero only in the case of soft hadron emission ||. 

Therefore the leading order expression for the cross section becomes 

K (i + (i-y) 2 ) 



dx dy dz 2y 2 



°i fi/p( x ) D h/i(z) 



i=q,q 



+ {l-x)M iMp {x,{l-x)z)\ . (3) 

In the next to leading order, the corrections to the four cross sections that can 
be defined taking into account the polarizations of the initial and final state 



1/2 spin hadrons, have been computed and can be found in refs. P, ID, III 



The scale dependence of fracture functions at 0(a s ) is driven by two kinds 
of processes, which contribute to the production of hadrons in the remnant 
target direction: the emission of collinear partons from those found in the 
target (the usual source of scale dependence of parton distributions, often 
called homogeneous evolution), and those where partons radiated from the 
one to be struck by the virtual probe, fragment into the measured hadron 
(the so-called inhomogeneous term). These two contributions lead at leading 
order to the following equation: 



dhgQ 2 v ' 2tt Jt/(i-0 u 



l m+ ° 77^ t?V> hi, M I Q' ) , (4) 




where Pj(u) and Pj' (u) are the regularized [|lj and real JRJ Altarelli-Parisi 
splitting functions, respectively. 

As was mentioned in the introduction, the ZEUS Collaboration has mea- 
sured DIS events identifying high-energy neutrons in the final state, at very 
small angles with respect to the proton direction {6i ab < 0.75 mrad), in the 
kinematical range given by 3 x 1CT 4 < x < 6 x 1CT 3 , 10 < Q 2 < 100 GeV 2 and 
high xl = E n / E p > 0.30. Here E p and E n are the target proton energy and 
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that of the produced hadron, but in the laboratory frame, which are straight- 
forwardly related to Eh and E p by the appropriate boost. The variables xl, 
9i ab and z are then correlated, specifically for 9 iab ~0,Xl~;z(1 — x). 

The ZEUS Collaboration have reported that events with xl > 0.50 repre- 
sent a substantial fraction (of the order of 10%) of DIS events, which means 
a contribution comparable in magnitude to the one given by the longitudinal 
structure function Fl- 

In the framework of a picture for semi-inclusive processes including frac- 
ture functions, as the one outlined above, the ZEUS findings can be repre- 
sented by 



Jo.50 dxdy P dx L dx L _ fo, 5 o M 2 " /P (x, X L , Q 2 ) dx L 



F$(x,Q 2 ' 



dx dy 



(5) 



an expression that is exact at leading order due to the fact that, as was 
mentioned earlier, at this order there is no contribution from currrent frag- 
mentation processes to the target region. In eq. (5) we have also defined the 
equivalent to F 2 for fracture functions: 

M 2 n/P (x, x L , Q 2 )=x^ e? Mi, n/P (x, x L , Q 2 ) (6) 

i 

and we have made explicit the integration over a finite (measured) range of 
Xl- Notice also that if the requirement 9\ ab ~ is not fulfilled, the exact 
relation between xl and z has to be taken into account. 

In fig. 1 we show the experimental outcome for this fracture function (as 
defined in eq. (5) and at Q 2 = 10 GeV 2 , taking advantage of the negligible 
Q 2 dependence of the data), and we compare it to F% and F[. We also show 
the contribution to the same observable coming from current fragmentation 
processes, computed taking into account the appropriate boost to the labo- 
ratory system and the cut in the 9i ab angle implemented in the experimental 
determination. It is worth noting that this last contribution is pure NLO and 
in fact it is about 8 orders of magnitude smaller than the experimental dataQ 
Furthermore, as can be seen in the figure, even releasing the constraint on 

1 For our analysis, we use distributions from refs. |l6|, [l7|] . 
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the angle for the final state neutron (but keeping the requirement of energy 
fraction > 0.50) the contribution still remains negligible, basically because 
of the kinematical suppression (the effect of the boost from the 7*p centre of 
mass system to the laboratory frame) and to the fact that this contribution 
begins at NLO. 

Fracture functions, as parton distributions in general, are essentially of 
a non-perturbative nature and have to be extracted from experiment. How- 
ever, their close relation with fragmentation and structure functions allows in 
certain extreme cases a model estimate for them. Basically, the task amounts 
to modelling the target fragmentation process and then writing down the re- 
sult in the language of fracture functions. This procedure provides input 
fracture functions, or at least hints of their functional shape, at a definite 
scale, which can then be compared with actual measurements and evolved 
with the corresponding evolution equations. 

Recently, a very sensible model estimate for the production of forward 
hadrons in DIS has been proposed [010, exploiting the idea of non-perturbative 



Fock components of the nucleon. In this approach the semi-inclusive DIS 
cross sections, and through them the corresponding fracture functions at a 
certain input scale Qq, can be interpreted as the product of a flux of neutrons 
in the proton (integrated over pj) times the structure function of the pion 
exchanged between them, i.e. 

M? /p (x, x L , Ql) ~ <f> n/p (x L ) Ft ( r ^-, Qg) . (7) 

Actually, ref. is a proposal to measure the structure function of pions 



at very small x using a non-perturbative computation of the flux, which is 
in very good agreement with experimental data on high energy neutron and 



A ++ production in hadron-hadron collisions |19 |. The structure functions 
of pions have been measured with some precision only for x > 0.10, but 
there are parametrizations based on dynamical parton model assumptions 
(valence-like distributions at low /Xg) ||20| , which have been shown an impres- 



sive predictive power at small x for the proton structure function |Tj|, and 



allow us to make a rough comparison between the model estimate and the 



2 In ref. |H( | an approach for the residual effects of fracture functions in the currrent 
fragmentation region has also been proposed. 
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data. In fig. 1 we show the that model estimate, taking Qq = 10 GeV 2 , 
is in excellent agreement with the data. In fig. 2 we also compare the xl 
dependence of the data and of the model prediction (normalized to the same 
number of events in the measured interval), showing also a good agreement. 
Regarding the shift in the position of the peak, it is worth noticing that the 
above-mentioned ZEUS data have not been corrected either for finite accep- 
tance effects or for those related with finite resolution, and the former may 
account for some, perhaps all, of the observed shift to largers values of xl 



21J| . This would eventually give an even better agreement between the the 



model and the data. 



3. Scale dependence of fracture functions: pions in the 
final state. 



The success of the model estimate encourages us to go further and use 
the functional dependence of fracture functions, induced by the model and 
corroborated by the data, to analyse also the Q 2 dependence. 

We first analyse the more familar process of neutron production. Since 
the probability of current parton fragmentation into a neutron (given by 
fragmentation functions) is comparatively small with respect to that of pro- 
cesses originated in the target (fracture functions), no significant effects are 
expected in the scale evolution arising from the inhomogeneous term in this 
process. The evolution is mainly driven by the usual homogeneous term 
of the evolution equations leading to an almost constant ratio between the 
number of neutron tagged events and that of all DIS events, as observed by 
ZEUS. 

However, the scale dependence induced in the cross section for the pro- 
duction of pions, at least in the kinematical region of very small x and small 
xl, can be considerably affected by the inhomogeneity, given that soft pions 
are produced more copiously from quarks than from neutrons^. Of course, 
in order to analyse these features of the evolution, an estimate for the pro- 
ton to pion fracture function at some input scale Qo ^ s required. For this 
purpose, we can use the same ideas formerly applied to neutron production, 

3 We are indebted to G. Veneziano for calling our attention to this point. 
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noticing two further advantages. First, the flux to be used can be straight- 
forwardly obtained from the one used in the last section by means of the 
crossing relation 

^+/ P (x L ) = ^/ p (l - x L ) , (8) 

where the superscript tt + indicates that the relation is valid only when the 
contribution due to the exchange of positively charged pions is considered. 
The second advantage is related to fact that the neutron structure function 
is much better known than that for pions. Taking into account these features 
and neglecting the insignificant contributions coming from A exchange, the 
proton to pion fracture function can be approximated by 

M? /p (x, x L , Ql) ~ ^ +/p (x L ) F? ( T -^— , Qg) . (9) 

Additionally, the direct measurement of this fracture function could be of 
interest as a double test of the crude factorization ideas. 

Notice that an analogous reasoning can be made considering the exchange 
and production of n° mesons, but with a proton instead of a neutron as 
final state and intermediate particle, respectively. The study of leading n° 
production may also be possible at ZEUS [f2lfl . 

In fig. 3 we show the model estimates for proton to pion fracture functions 
(taking Ql = 4 GeV 2 ), integrated over two different bins of xl, compared 
with the contribution coming from the current fragmentation processes. We 
assume here the same restrictions as in the data from the ZEUS Collabo- 
ration for neutron production^. As can be seen, higher order fragmentation 
contributions are quite less significant than fracture contributions, even at 
the lowest xl bin. 

Of course, the model is not expected to work over the whole kinematical 
region and, in fact, any deviation from the scale dependence implied in eq. 
(9) (note that the flux is assumed to be Q 2 independent) would show the 
breakdown of the approximated factorization hypotesis. However, the ansatz 
in eq. (9) can be taken as an effective relation, valid at some initial value 

4 Actually, it would also be interesting to count with the pt distribution of for- 
ward hadrons in order to establish the exact range where fracture function contribution 
dominates. 
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of Qq, for which the estimated flux is adequate, and therefore provides a 
sensible input distribution^. As usual, the correct scale dependence is that 
given by the evolution equations for fracture functions, and that is the aim 
of our next step. 

In order to study the effect of the inhomogeneity in the evolution we take 
different values of xl, and keep them fixed while we analyse the x and Q 2 
dependence of fracture functions induced by both the homogeneous and the 
complete evolution equations. 

In fig. 4 we show the result of an evolution from Qq = 1 GeV 2 at 
xi — 0.50. Both solutions (the homogeneous and the complete) are super- 
imposed, the difference being less than 0.1%. These behaviours are perfectly 
compatible with the results obtained by the ZEUS Collaboration in the case 
of neutron production (where the inhomogeneity contributes about 10 times 
less) in the same kinematical region, where no difference has been found in 
the evolution between F2 and M2. 

However, for smaller xl the situation is completely different. As the 
fragmentation function increases with lower values of the argument, the in- 
homogeneous contribution becomes much more relevant and its effect in the 
evolution is sizeable. In fact, fig. 5 shows the evolution result for xl = 0.10 
and Ql = 1 GeV 2 , where the full evolution results outsize the homogeneous 
one by a factor of 4 at small x. These corrections are smaller if the ansatz of 
eq. (9) is assumed to be valid at values of Ql = 4 GeV 2 (fig. 6) and Ql = 10 
GeV 2 (fig. 7), but still remain considerable and show the same behaviour. 

The uncertainty due to the estimate of the input fracture functions at 
some initial scale obviously will be lifted as soon as they are measured; it 
will be then possible to study their evolution in a model-independent way, 
completely determined by QCD, and to compare the result with the experi- 
mental data at different values of the scale. 

In fig. 8, we show for completeness the results obtained using a bin in 
Xl between 0.10 and 0.20. In that case the weight of the corrections is low- 
ered by the inclusion of larger-x^ terms, but it is noticeable that they are 
still about 20%, whereas corrections from the fragmentation processes to the 

5 In fact, different choices of Qq do not modify the overall trends of the evolution. 
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fracture function, as in fig. 3, are less than 1% in the small-x range. This 
is due to the kinematical restrictions imposed over the final state, but also 
because to the fact that the fragmentation contribution appears at the same 
order of a s in the evolution, but one order higher in the cross section than 
that coming from the fracture one. 

4. Summary and conclusions 

In this paper, we have analysed recent experimental data on the produc- 
tion of forward neutrons in DIS in terms of fracture functions, finding that 
the main features of the data can be fairly reproduced by this perturbative 
QCD approach, once a non-perturbative model estimate for the input frac- 
ture functions is given. Studying the evolution properties of these fracture 
functions in the specific case of forward pions in the final state, we have 
found that the effects of the inhomogeneous term in the evolution equations 
are large and measurable, particularly in the kinematical region of very small 
x (accessible to HERA) and small xl- These effects are negligible for large 
values of xl, justifying the use of the usual homogeneous Altarelli-Parisi 
equations for, as an example, the t-integrated diffractive structure function 
F^ 3 \x p} (3, Q 2 ), which is just the fracture function of protons in protons 
M% /P ((3x p , (l-x p ),Q 2 ) §i. 

The measurement of the production of pions in the target fragmentation 
region in DIS will be useful to check the range of applicability of factorization 
ansatze for fragmentation processes and to perform a new test of QCD by 
means of the analysis of the evolution of fracture functions, which is predicted 
to be different from the usual one in the kinematical region alluded to before. 
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Figure 1: Fracture function of neutrons in protons as measured by ZEUS 
compared with the model prediction and current fragmentation contributions 



100 




0.2 0.4 0.6 0.8 1 1.2 

Figure 2: The xl dependence of the model compared with the ZEUS data, 
integrated in the region 3 x 10 -4 < x < 6 x 
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Figure 3: Prediction for the fracture function of ir + in protons for two differ- 
ent bins of Xl and the current fragmentation contribution 




15 




16 




17 



